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S  WARY 


It  is  shown  that  non-negatlvlty  of  the  solution  of  the 
heat  equation,  given  non-negative  Initial  values,  and  suitable 
boundary  conditions,  can  be  established  quite  readily  once  the 
existence  of  a  solution  of  the  equation  depending  continuously 
upon  the  Initial  values  has  been  demonstrated. 

It  Is  shown  that  this  property  Is  trivially  true  for  the 
solution  of  the  appropriate  finite  difference  approximation  to 
the  partial  differential  equation,  and  that  the  convergence  of 
the  solution  of  the  finite  difference  equation  to  the  solution 
of  the  original  equation  Is  quite  simple,  under  the  above 
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ON  THE  NON-lfEOATIVITY  OP 
SOLUTIONS  OP  THE  HEAT  EQUATION 

Richard  Bellman 


t^I .  INTRODUCTION 

In  treating  a  functional  equation,  once  the  question  of 
existence  and  uniqueness  has  b^en  disposed  of,  we  turn  to  a 
more  precise  study  of  the  analytic  character  of  the  solution. 

It  frequently  happens  that  a  method  which  works  very  efficiently 
to  establish  existence  and  uniqueness  does  not  yield  other 
properties  of  the  solution  In  any  ready  fashion.  Conversely, 
methods  which  yield  non— negativity ,  convexity,  and  so  forth, 
may  not  be  Ideally  suited  for  the  establishment  of  the  basic 
properties.  However,  a  combination  of  several  techniques  may 
yield  the  results  we  desire  quite  easily. 

To  Illustrate  these  remarks,  let  us  consider  the  heat 
equation 


“t  ■  “xx 

(l)  u(x,0)  -  v(x),  0  <  X  <  1, 

u(0,t)  ■  u(l,t)  -  0,  t  >  0. 

We  shall  assume  that  we  have  demonstrated,  by  some  means  or 
other,  the  existence  of  a  solution  which  depends  continuously 
upon  v(x)  In  the  L  — norm  for  t  >  0,  and  that  this 
solution  Is  unique.  As  we  shall  see,  a  method  based  upon 
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flnlte  differences  will  enable  us  to  demonstrate  the  fact  that 
this  solution  Is  non-negative  for  t  >  0,  provided  that  It  Is 
non-negative  at  t  -  0,  l.e.  provided  that  v(x)  >  0,  for 
0  <  X  <  1.  On  the  other  hand,  the  particular  proof  used  to 
establish  existence  and  uniqueness  may  not  have  yielded  non— 
negativity  In  a  simple  fashion,  and,  as  Is  known,  an  existence 
and  uniqueness  proof  based  upon  finite  differences  Ic  not  a 
completely  simple  matter. 

“t  ■  “xx 

To  Illustrate  our  Ideas,  begin  with  the  simpler  equation 
“t  ■  “xx 

amd  consider  the  difference  scheme 

(a)  w(x,t  +  5^2)  -  [w(x-»^,t)  -»■  w(x-^,t)]/2, 

(2) 

(b)  w(x,0)  -  v(x), 

where  x  takes  the  values  0,  6,  26,  ...,  1,  and  t  assumes 

2  2 

the  values  0,  6  /2,  6  ,  ...  .  The  function  w(x,t)  Is 
defined  by  linearity  at  non— lattice  points. 

It  Is  easy  to  see  that,  formally,  the  recurrence  relation 
approaches  the  partial  differential  equation  as  6  — >  0. 

As  mentioned  above,  a  rigorous  proof  that  the  solution  of 
(2)  converges  to  the  solution  of  (l)  as  6  — >  0,  starting 
from  first  principles.  Is  non— trivial.  However,  as  we  shall 
see  below,  a  proof  of  this  fact  Is  quite  simple,  once  we  have 
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established  the  existence  and  uniqueness  of  a  solution.  The 
fact  that  w(x,t)  Is  non-negative  for  any  8  >  0  Is  Immediate, 
and  this  yields  the  conclusion  that  u(x,t)  >  0. 

Since  we  have  assumed  the  existence  of  a  solution  of  (l) 
which  Is  a  continuous  function  of  v(x),  there  is  no  loss  of 
generality  In  assuming,  for  our  current  purposes,  that  v(x) 
possesses  appropriate  continuity  properties,  sufficient  to 
ensure  that 


(3) 


'^xxxx^^'^^ 


<  m  <  00  , 


where  R  Is  the  bounded  region  0<t<T<(r-. 

We  may  for  example  take  v(x)  to  be  a  trigonometric  poly¬ 
nomial.  Under  the  assumption  of  (2),  It  Is  easy  to  show  that 
In  R 

(^)  11m  w(x,t)  -  u(x,t ) . 

6—^  0 

We  have,  by  virtue  of  (2) 

(5)  u(x,t  8^/2)  ■  [u(x'»*6,t)  u(x— 6,t)]/2  +  6^r(x,t), 

where  |r(x,t)|  <  2in  in  R.  Consequently,  the  function 
z(x,t)  ■  w(x,t)  —  u(x,t)  satisfies  the  recurrence  relation  of 
(la)  with  the  Initial  condition  z(x,t)  -  0.  Let 

(6)  d(t)  -  Max  |w(x,t)  -  u(x,t)|. 

0<x<l 
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Using  the  recurrence  relation,  we  see  that 

(7)  d(t-M5V2)  <  d(t)  +  26^m 

p  2 

for  t  -  0,  6  /2,  6  ,  whence 

(8)  d(t)  <  2fi^mN 

for  0  <  t  <  6^N.  Let  6^N  -  T.  Then 

(9)  d(t)  <  26^mT 

for  X  and  t  In  R. 

Prom  this,  we  obtain  the  desired  result  as  6  — >  0. 

Note  that  we  can  also  conclude  from  the  foregoing  result 
that  u(x,t)  Is  concave  In  x  for  any  value  of  t  If  v(x) 
Is  concave  In  x.  This,  In  turn.  Implies  that  u(x,t)  Is 
decreasing  In  t  for  each  fixed  value  of  x. 

'^t  ■  '^xx 

To  extend  the  same  argument  to  the  general  equation  of 
(l.l),  we  employ  the  recurrence  relation 

w(x,t  +  6^/2)  -  + 

2 

X  -f  q(x,t  )ft^/2  /  4.  \  . 

,  [  o  u(y,t)dy. 

-  q(x,t)6^/2 


(1) 
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If  W'-'  assume  that  q(x,t)  >  0  for  0<x<l,  t>0, 

the  recurrence  relation  above  shows  that  w(x,t)  ^  ^ 

X  and  t.  The  proof  that  w(x,t)  converges  to  u(x,t)  as 
0  — >0  follows  the  same  .  tnes  as  befc  ♦* . 

To  see  that  It  Is  sufficient  to  assumt;  ‘'at 


(2)  q(x,t)>-X>-oo,  0<x<l,  0^’  T, 

for  any  T  >  0,  where  X  ■  ^'(T),  we  proceed  as  follows. 
Write 

(3)  u  -  e”^^v. 

Then  the  equation  of  (l.l)  becomes 

(4)  v^  -  v^^  +  (q(x,t )-hX)v, 

with  the  same  boundary  conditions.  The  new  function 

(5)  qi(x,t)  -  q(x,t)  -h  X 


Is  non-negative. 

.  GENERALIZATIONS 

It  Is  clear  that  the  same  method  may  be  employed  to  obtain 
corresponding  non-negatlvlty  results  for  the  solution  of  the 
heat  equation  for  higher  dimensions  smd  arbitrary  regions.  The 
essential  part  of  the  proof  Is  the  a  priori  demonstration  of  the 
existence  and  uniqueness  of  a  solu' ion  depending  continuously 
upon  the  Initial  values,  In  an  apj roprlate  metric. 

Similarly,  a  number  of  corri'spondlng  resulte  ''an  be 
established  for  various  classes  of  nonlinear  equations. 
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